Evanescent waves in photonic crystals and image of Veselago lens 
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It is shown that negative electric permittivity e and magnetic permeability fi recently discovered 
in a photonic crystal in the vicinity of the F-point are properties of propagating modes only. The 
evanescent modes rather decay than increase in the bulk of the crystal though they may be amplified 
by surface waves. If surface support such waves, the evanescent waves may improve the image of a 
thin Veselago lens. It is shown that a "perfect lens" contradicts to the wave optics and a criterion 
of "superlensing" is formulated. 
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Left handed medium (LHM), defined by Veselago pj as 
a medium with simultaneously negative \i and e, recently 
attracted much attention mostly because of the negative 
refraction at the interface with a regular medium (RM). 
This effect allows creation of a unique device called the 
"Veselago lens". This lens is a slab of LHM inside a 
RM with a condition that both media have the same 
isotropic refractive index and the same impedance. The 
interest in LHM was significantly amplified by the work 
of PendryQ who argued that the Veselago lens is a "per- 
fect lens" in the sense that it gives a perfect image of 
the point source. This statement is based upon the ob- 
servation that the evanescent waves (EW's) of a form 
exp(ik y y — kx) that usually decay in the near field re- 
gion are amplified by the LHM. Pendry claimed that the 
amplified EW's may restore a perfect image. Follow- 
ing Veselago, Pendry considered a "hypothetical" LHM 
(HLHM) with negative and real fj, and e. He did not 
present a solution in a coordinate space. This was done 
by Ziolkowski and HeymanQ and their solution reveals 
a fundamental problem with Pendry's idea. The solu- 
tion diverges exponentially at each point of a 3D domain 
near the focus, just where fields of EW's increase due to 
amplification. This was found out almost simultaneously 
by three groups 0,0,0- The first group argued that the 
EW's should just be omitted as non-physical. The sec- 
ond group used diffraction theory, known to be exact at 
small wave length, and found a finite width of the focus. 
Haldanejgt argued that the problem does not have a so- 
lution in the framework of macroscopic electrodynamics. 

We consider here a dielectric photonic crystal (PC) 
that is known to be a real LHM in some frequency 
range0 H}. We show that its electrodynamics differs 
from that of the HLHM and that in the real system neg- 
ative fi and e of propagating waves do not provide am- 
plification of EW's. We present also a general proof that 
the perfect imaging of a point source contradicts wave 
optics independent of the nature of a lens. 

A large group of works pioneered by Luo et al. 0, 0] 
claims observation of negative refraction, focusing, and 
even superlensing in a PC due to photons with wave vec- 
tors deep in the Brillouin zone. Propagation of such pho- 
tons can be described macroscopically only taking into 
account spatial dispersion (k-dependence of e) . The the- 



ory of Veselago can hardly be generalized for the case of 
spatial dispersion mostly because of an extra term in the 
definition of Poynting's vector [ll) . Following Ref.0 we 
do not consider such a medium left handed and do not 
discuss it here. 

It was shown recently that a 2D dielectric uniax- 
ial PC made of non-magnetic materials can behave as a 
LHM with negative e and /i if it has a negative group ve- 
locity in the vicinity of the T-point. This was proved for 
propagating modes only. Experimental demonstration of 
negative refraction in a metallic PC usin g th e modes near 
the T-point has been presented in Ref. |l2j. 

In this paper we consider p-polarization for both EW's 
and propagating modes in a uniaxial PC. Fig. |T^a) shows 
the elementary cell in the plane x — y and the spectrum 
of p-polarized propagating waves. The dashed line is our 
working frequency. All data below are obtained for this 
frequency. 

From the same numerical solution of microscopic 
Maxwell's equation that gives the spectrum one gets mi- 
croscopic fields h z (x, y), e y {x, y) and e x {x, y) correspond- 
ing to a given frequency of propagating modes. The mi- 
croscopic fields are the Bloch functions wile the macro- 
scopic fields have a form B z —<h z >, E =<e>, where 
<...> means averaging over the unit cell. Near the T- 
point the Bloch functions have a small k, and the macro- 
scopic fields of the propagating modes are plane waves 
that obey the macroscopic Maxwell equations. 

As follows from Ref.[8j, the values of \l zz = [i and 
txx — £yy = Ej that describe propagation of these waves, 
are negative because their group velocity is negative. 
Here z is the axis of the crystal. It has been shown 
previously0 that the values of e and fj, can be expressed 
through the integrals of Bloch's function at k — 0. It is 
important that in a system with a spatial dispersion the 
separation of magnetization and displacement currents is 
impossible Our parameters e(u>) and (x(oj) are not 

exclusive properties of the medium at a given frequency, 
as in macroscopic electrodynamics without spatial dis- 
persion. Rather they are properties of a particular mode 
of this medium, whose fields are used for the calcula- 
tions. We argue below that they are not applicable to 
the EW's. The theorem^ connecting signs of both e(w) 
and fi(uj) with the sign of group velocity is not applicable 
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to the EW's as well. Near the T-point ~ & 2 0. Since 
the dispersion law u>(k) is also isotropic one obtains fj, 
that are function of ix> only. This approach is good for 
an infinite PC, but may create problems with boundary 
conditions [ll| because in coordinate space k 2 should be 
considered as Laplacian operator. We are going to re- 
fer to this problem somewhere else. For the EW's this 
approach fails completely because their dispersion law is 
anisotropic. 

To demonstrate Veselago lens one should know fi and e 
of the RM medium around the PC. In this paper we find 
the value of of refractive index n — y/ep from the function 
u>(k), shown in Fig. [Ja), using definition u> — kc/n. The 
second condition that provides matching of impedance is 
the absence of reflection of the incident plane wave in a 
wide range of incident angles. 

Using these two conditions we get that at the working 
frequency the RM material should have e = 1.125, fi — 
0.08. The physical reasons for the appearance of /i in a 
non- magnetic PC are discussed in Ref. Q. In all com- 
putations below we use these values for homogeneous re- 
gions surrounding the PC slab. 

Fig. ^b) shows the result of our simulation for light 
propagation through the PC slab surrounded by a homo- 
geneous medium. One can see the negative refraction of 
light coming in and going out of the PC with equal abso- 
lute values of incident and refracted angles and without 
any visible reflections from the interfaces. We checked it 
in a wide range of the incident angles (from to 60 de- 
grees). An animation^! shows that the wave fronts are 
moving to the right outside the PC slab and to the left in- 
side the slab. All of these phenomena exactly correspond 
to Veselago 's theory Q. 

Now we come to EW's in the same PC slab imbed- 
ded into the same RM as in Fig 1 (b) . Our results 
are shown in Fig. 2. The boundary conditions implied 
for an incident wave at x = are h z = exp (ik y y) and 
dhz/dx = —Kh z - Note that due to reflection the total 
field shown in Fig. 2 may not coincide at x — with 
the incident one. To obey Maxwell's equation we imply 
the dispersion law of EW's in a form k 2 — k 2 — uj 2 n 2 /c 2 . 
The distance from the left boundary of the RM (point 
x = 0) to the PC slab is 1/2 of the slab thickness. In 
Fig. 2 (a)the PC slab is substituted by the HLHM that 
has negative e and fj, with the absolute values same as in 
the RM. One can see a strong amplification of the EW 
first discovered by Pendry 2]. However the PC slab gives 
a very different picture. The result strongly depend on 
the way the surface is cut. Our surfaces are perpendicu- 
lar to [10] direction, but they cut either across the holes 
(AH) or between the holes (BH). In the case of the BH 
surface (dashed lines) the EW decreases in the PC slab 
in the same way as in RM without any signature of am- 
plification. However, there is an amplification at the left 
AH surface for a thick slab and at both AH surfaces for a 
thin slab. We think that this amplification are due to the 
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FIG. 1: (Color online) a) Five lowest bands of the photonic 
spectrum of a 2D PC with a period d. The unit cell of the PC 
is shown in the inset. The PC consists of a square lattice of 
circular cylindrical holes in a dielectric matrix with e m = 12, 
|U m = 1. The radii of the holes R = 0.35d. b) Negative re- 
fraction of light propagating through the PC slab surrounding 
by homogeneous medium at the working frequency. Periodic 
boundary conditions are used in the vertical direction. The 
arrows show the directions of the wave vector. 



surface waves (SW), excited by EW's. Since the picture 
is not symmetric with respect to left and right interfaces 
this excitation is not resonant. These SWs have noth- 
ing in common with left-handed properties and with po- 
larons in the HLHM, considered by Ruppin^] because 
they depend on the properties of the surface. We think 
these SWs appear just due to the termination of peri- 
odic structure. They have been studied before by many 
authors[l3 H [H llSW with and without connection to 
the focusing. 

Thus, we have found that EW's cannot be described 
in terms of the same negative e and [i that are found 
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FIG. 2: (Color online) Magnetic field of the EW in the 
HLHM (a) and in the PC slabs (b),(c) surrounded by a regu- 
lar homogeneous medium at the following parameters of EW: 
k y — (\/5/2)ko, k — 0.5fc , fco = tun/c, n = ^fejl = 0.3. The 
cross section y — is shown. The thickness of the PC slab is 
20d (b) and lOd (c). The full lines (black) are for the AH sur- 
faces, the dashed lines (red) are for BH surfaces. The period 
of oscillation is d. 



for the propagating waves. The formal reason is that 
the dispersion law of EW is anisotropic even near the V- 
point. Therefore in this case fj, cannot be represented in 
k-independent form as in the case of propagating waves. 
The physical explanation reads that the EW in the direc- 
tion of decay is the same for both regular and backward 
wave. Therefore in this direction there is no difference 
between LHM and RM. We have verified, however, that 
if an EW propagates (has a real fc) in a direction not 
parallel to the surface of the PC slab, it shows a negative 
refraction, similar to that shown in Fig. 1(b). 
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FIG. 3: (Color online) Distribution of magnetic energy near 
the focus of the lens in lateral and perpendicular directions. 
Solid lines show results of simulation. The upper curve (red) 
is distribution in lateral direction for a thin lens with AH 
surfaces, the lower curve(blue) is the same for a thick lens. 
Open circles (red) shows the analytical result for a thin lens, 
the dashed curve (blue) shows the same for a thick lens. Dash- 
dotted line (black) shows the distribution of magnetic energy 
of the point source. The inset shows computational distri- 
bution in perpendicular direction (solid blue) and analytical 
result in this direction (dashed blue) both for the thick lens 
only. 



Now we describe the image of the Veselago lens. 
The point source is the 2D Green function H s = 
iHoHq (pko) exp —icot, where fco = tun/c and the Hankel 
function Hq — Jq + iNq. It is located in the RM at a 
distance a to the left of the PC slab. 

It is easy to show that that the perfect image of this 
source is impossible. Indeed, if the image has the same 
form as the source, the magnetic field near the focus 
would obey inhomogeneous Helmholtz equation with a 
delta-function at the focus, while the focal point would 
not contain any source. Thus, this solution does not obey 
homogeneous Helmholtz equation near the focal point. 
Note that this argument docs not specify the nature of a 
lens. 

We simulated the Veselago lens with the same PC slab 
and the same parameters of the RM. Similar to geom- 
etry of Fig. 2 the distance a is 1/2 of the thickness of 
the PC slab. Fig. 3 shows the spatial distributions of 
magnetic energy near the focus in both lateral (y) and 
perpendicular (x) directions. The distribution in the lat- 
eral direction is symmetric for y and —y so it is shown 
for positive y only. 

Our results shown in Fig. 2 give a possibility to cal- 
culate analytically the distribution of energy near the 
focus of a thick Veselago lens with both AH and BH sur- 
faces or of a thin lens with BH surface. In these cases 
evanescent waves do not reach the focus. At positive x 
the Green function of the source can be represented in a 
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form H s = H v + H ev , where 
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contains only EW's. 

If the EW's decay inside the PC slab, only term H p 
contributes to the focus if the slab is thick enough . On 
the other hand our results show that for the propagat- 
ing modes the PC slab works exactly like a HLHM. Us- 
ing the Fresnel equation it is easy to showQ that the 
field near the focus is given by H p (x',y'), where the ar- 
guments originate at the focal point. This function is 
not symmetric with respect to x' and y': the width of 
the distribution in perpendicular to the slab direction is 
larger than in the lateral direction (See Fig. This gen- 
eral result originates from the withdrawal of EW's. Note 
that the source field depends only on p. Another im- 
portant feature of this imaging is that the image in the 
lateral direction is H p (0,y') = iHoJo(y'ko)exp—icjt. If 
we consider real fields, the part of the source field that 
is proportional to sin(u;t), namely H J (pko) sin(wi), has 
a perfect image in the lateral direction while the part of 
the source with cos(wi) has a zero image in this direction. 
This follows only from the plane (non-axial) geometry of 
the lens that is assumed to be infinite in the lateral direc- 
tion. The perfect lateral imaging of the Bessel function 
Jo does not contradict to general principles because this 
Bessel function obeys homogeneous Helmholtz equation. 

In our simulation, however, the slab has a finite length 
h in y-direction. In this case some propagating modes 
go outside the slab. To take this into account one 
should integrate in Eq. Q from — k m to k m , where 
k m = ko(h/2a)/ yT + (h/2a) 2 , a is the distance from the 
source to the slab. In our case h/2a = 4, 8 for thick and 
thin lenses respectively. Due to this factor the lateral 
distribution slightly differs from Jq. These calculations 
are shown in Fig. 3 for both thick (20d) and thin (lOd) 
lenses. They are different because of different values of 
a. For the thick lens analytical calculation is in a very 
good agreement with the computational data. 

A sharp first peak of a Bessel-like function in 
lateral direction is often considered as a result of 
"superlensing" [2(| • We think, however, that any re- 
sult that follows from Eq.Q has nothing to do with "su- 
perlensing" . Following Pendry we would define super- 
lensing as a result of amplification of EW's, but Eq.0J 
does not contain EW's at all. It gives just a diffraction 



pattern of the Veselago lens that was obtained in Ref. 5] 
for a 3D-lens. Then the energy distribution given by this 
equation should be called "regular lensing" . Thus, our 
thick Veselago lens based upon PC does not provide any 
superlensing. The thin lens with AH sufrace gives sub- 
stantially sharper image than it follows from Eq. with 
a proper /i/2a-correction(See Fig. 3). This is due to am- 
plification of EW's shown in Fig. 2(c) and this effect can 
be considered as a superlensing. 

In conclusion we have shown that the Veselago lens 
imaging has some peculiar features, such as a perfect 
imaging of the phase shifted part of the source. We show 
that a thin lens in a near field regime may provide a 
substantial superlensing depending on the way the slab 
is terminated. The physics of this superlensing is not 
connected with the left-handed properties of a medium. 

The work has been funded by the NSF grant DMR- 
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